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ABSTRACT 


Context. To understand the origin and formation pathway of wide-orbit gas giant planets, it is necessary to expand the limited sample 
of these objects. The mass of exoplanets derived with spectrophotometry, however, varies strongly as a function of the age of the 
system and the mass of the primary star. 

Aims. By selecting stars with similar ages and masses, the Young Suns Exoplanet Survey (YSES) aims to detect and characterize 
planetary-mass companions to solar-type host stars in the Scorpius-Centaurus association. 

Methods. Our survey is carried out with VLT/SPHERE with short exposure sequences on the order of 5 min per star per filter. The 
subtraction of the stellar point spread function (PSF) is based on reference star differential imaging (RDI) using the other targets 
(with similar colors and magnitudes) in the survey in combination with principal component analysis. Two astrometric epochs that 
are separated by more than one year are used to confirm co-moving companions by proper motion analysis. 

Results. We report the discovery of YSES 2b, a co-moving, planetary-mass companion to the K1 star YSES 2 (TYC 8984-2245-1, 
2MASS J11275535-6626046). The primary has a Gaia EDR3 distance of 110 pc, and we derive a revised mass of 1.1 Mo and an age 
of approximately 14 Myr. We detect the companion in two observing epochs southwest of the star at a position angle of 205? and 
with a separation of ~ 1705, which translates to a minimum physical separation of 115 au at the distance of the system. Photometric 
measurements in the H and K; bands are indicative of a late L spectral type, similar to the innermost planets around HR 8799. We 
derive a photometric planet mass of 6.3756 Myyy using AMES-COND and AMES-dusty evolutionary models; this mass corresponds 
to a mass ratio of q = (0.5 + 0.1) % with the primary. This is the lowest mass ratio of a direct imaging planet around a solar-type star 
to date. We discuss potential formation mechanisms and find that the current position of the planet is compatible with formation by 
disk gravitational instability, but its mass is lower than expected from numerical simulations. Formation via core accretion must have 
occurred closer to the star, yet we do not find evidence that supports the required outward migration, such as via scattering off another 
undiscovered companion in the system. We can exclude additional companions with masses greater than 13 Mj, in the full field of 
view of the detector (0715 < p < 5750), at 075 we can rule out further objects that are more massive than 6 Mj, , and for projected 
separations p > 2” we are sensitive to planets with masses as low as 2 Mjyp . 

Conclusions. YSES 2b is an ideal target for follow-up observations to further the understanding of the physical and chemical forma- 
tion mechanisms of wide-orbit Jovian planets. The YSES strategy of short snapshot observations (< 5 min) and PSF subtraction based 
on a large reference library proves to be extremely efficient and should be considered for future direct imaging surveys. 


Key words. Planets and satellites: detection — Planets and satellites: formation — Instrumentation: high angular resolution — 
Techniques: image processing — Stars: individual: YSES 2 (TYC 8984-2245-1) 


1. Introduction 


Despite several remarkable exoplanet and brown dwarf discov- 
eries by high-contrast imaging at high angular resolution in the 


* Based on observations collected at the European Organisa- 
tion for Astronomical Research in the Southern Hemisphere under 
ESO programs 099.C-0698(A), 0101.C-0153(A), 0101.C-0341(A), and 
106.20X2.001. 


past few years (e.g., Marois et al. 2008; Schmidt et al. 2008; 
Marois et al. 2010; Lagrange et al. 2010; Rameau et al. 2013; 
Bailey et al. 2014; Macintosh et al. 2015; Chauvin et al. 2017; 
Keppler et al. 2018; Haffert et al. 2019; Janson et al. 2019; Bohn 
et al. 2020a,b), there is an ongoing debate regarding the forma- 
tion mechanisms that create these super-Jovian gas giants with 
semimajor axes greater than 10au. It is unclear whether these 
companions have a star-like origin from a collapsing molecu- 
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lar cloud that is broken up into fragments, creating planetary- 
mass objects similar to a stellar binary (Kroupa 2001; Chabrier 
2003), or through formation in a circumstellar disk instead. The 
classical bottom-up framework postulates formation via core ac- 
cretion by coagulation of small dust grains into planetary em- 
bryos (Pollack et al. 1996; Alibert et al. 2005; Dodson-Robinson 
et al. 2009; Lambrechts & Johansen 2012). These evolve either 
via collisions or pebble accretion (Johansen & Lacerda 2010; 
Ormel & Klahr 2010) into planetary cores that are massive 
enough to accrete a gaseous envelope and to open a gap in the 
disk (Paardekooper & Mellema 2004). In the corresponding top- 
down scenario, planetary cores can be created by gravitational 
instabilities leading to the collapse of dense regions in the pro- 
toplanetary disk (Boss 1997; Rafikov 2005; Durisen et al. 2007; 
Kratter et al. 2010; Boss 2011; Kratter & Lodato 2016). 


To study this question for the underlying planet forma- 
tion mechanisms from a statistical point of view, several direct 
imaging surveys have been conducted (e.g., Vigan et al. 2012; 
Galicher et al. 2016; Bowler 2016; Vigan et al. 2017). Synthetic 
planet populations that represent each of the potential formation 
channels can be compared to the observational results from sur- 
veys and place constraints on the efficiency of the correspond- 
ing formation pathway (e.g., Mordasini et al. 2009a,b; Forgan & 
Rice 2013; Forgan et al. 2018). The two largest surveys were car- 
ried out with two of the most advanced adaptive-optics assisted, 
high-contrast imagers available: the Spectro-Polarimetric High- 
contrast Exoplanet REsearch (SPHERE; Beuzit et al. 2019) in- 
strument at the 8.2m ESO/VLT and the Gemini Planet Im- 
ager (GPI; Macintosh et al. 2014) at the 8.1 m Gemini South 
telescope. The preliminary statistical analysis of the first 300 
stars from the Gemini PLanet Imager Exoplanet Survey (GPIES; 
Nielsen et al. 2019) concludes that giant planets between 10 au 
and 100 au that have masses smaller than 13 Mj, favorably form 
via core accretion mechanisms, whereas brown dwarf compan- 
ions in the same separation range but with masses from 13 Mjyp 
to 80 Mjyp seem to be predominantly created by disk instabil- 
ities. This finding is supported by the analysis of the first 150 
stars observed within the scope of the SpHere INfrared survey 
for Exoplanets (SHINE; Vigan et al. 2020), which additionally 
hypothesizes that companions with masses between 1 Mj, and 
75 My, are likely to originate from bottom-up formation sce- 
narios around B and A type stars, whilst objects of the same 
mass around M-type stars are consistent with simulated popula- 
tions from top-down mechanisms. For the intermediate masses 
of F-, G-, and K-type stars, the observed detections can be ex- 
plained by a combination of both formalisms. A statistical meta- 
analysis on the distribution of wide-orbit companion eccentrici- 
ties carried out by Bowler et al. (2020) provides supporting ev- 
idence for two distinct formation channels shaping the popula- 
tions of giant planets (2 Mjup < M < 15 Mjyp) and brown dwarfs 
(15 Mygp < M < 75 Miop). 


Most of these statistical evaluations are affected by the small 
number of actual substellar companions that were detected in 
the preceding imaging surveys. To expand the sample size for 
solar-mass host stars, we started the Young Suns Exoplanet Sur- 
vey (YSES; Bohn et al. 2020a, Bohn et al. in prep.) that is ob- 
serving a homogeneous sample of 70 ~15 Myr-old, K-type stars 
in the Lower Centaurus Crux (LCC) subgroup of the Scorpius- 
Centaurus association (Sco-Cen; de Zeeuw et al. 1999). All stars 
have masses close to 1 Mo and the proximity (average parallac- 
tic distance (D) = 114 + 17; Gaia Collaboration et al. 2020) and 
youth of the LCC facilitate the direct imaging search of giant, 
self-luminous substellar companions around these stars. 
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In this article we report the detection of a new exoplanet that 
was discovered within the scope of our survey. As this is already 
the second planetary system discovered by YSES, we introduce 
a new stellar identifier that is based on our survey acronym. The 
details of this new designation are described in $2. In $3 we 
describe our observations and data reduction methods. We dis- 
cuss previous observations on the host star and reassess its main 
parameters in $4. The results of our high-contrast imaging ob- 
servations are presented and analyzed in $5. We discuss poten- 
tial formation mechanisms of this newly detected exoplanet in 
$6 and we present our conclusions in $7. 


2. Nomenclature of YSES planets 


Owing to the recent success of YSES, we decided to introduce 
a dedicated catalog that will be used for star-planet systems dis- 
covered within the scope of our survey. The YSES acronym has 
been verified by the IAU Commission B2 Working Group on 
Designations and was added to the Simbad database (Wenger 
et al. 2000).! The nomenclature of planet hosts from our survey 
is YSES NNN and planets that are associated with these stars 
will be named YSES NNNa, accordingly. Following these guide- 
lines, we assigned the host star of the intriguing multi-planet 
system that was discovered around TYC 8998-760-1 the new 
primary identifier YSES 1 (Bohn et al. 2020a,b). The planets 
formerly known as TYC 8998-760-1 b and TYC 8998-760-1 c, 
will be named YSES 1b and YSES 1c, henceforth. Further plan- 
etary systems discovered by our survey will receive designated 
YSES identifiers followed by ascending integer identifications 
(IDs). Hence, the new companion discovered within the scope 
of this paper will be referred to as YSES 2b, orbiting its Sun-like 
host YSES 2. 


3. Observations and data reduction 


We observed YSES 2 (TYC 8984-2245-1, 2MASS J11275535- 
6626046) as part of YSES on the nights of 2018 April 30 (PI: 
Kenworthy) and 2020 December 8 (PI: Vogt) with SPHERE 
(mounted at the Nasmyth platform of Unit Telescope 3 of the 
ESO Very Large Telescope). We used the IRDIS camera (Dohlen 
et al. 2008) in classical imaging mode, applying a broadband fil- 
ter in the H and K; bands during the first and second nights, re- 
spectively. The observations were carried out in pupil stabilized 
imaging mode and an apodized Lyot coronagraph was used to 
block the flux of the primary star (Soummer 2005; Martinez et al. 
2009; Carbillet et al. 2011). In addition to the science frames, we 
obtained center frames, with a sinusoidal pattern applied to the 
deformable mirror that creates a waffle pattern to locate the posi- 
tion of the star behind the coronagraphic mask; sky frames of an 
offset position with no adaptive optics (AO) correction and with- 
out any source in the field of view, to subtract the instrument and 
thermal background; and non-coronagraphic images of the star 
that are used for photometric reference of point sources detected 
in the science images. This last category of non-coronagraphic 
flux images was obtained with an additional neutral density filter 
in the optical path to record an unsaturated stellar point spread 
function (PSF) in the linear readout regime of the detector. This 
neutral density filter (filter ID: ND 1.0) provided an attenuation 
of 7.9 and 6.9 across the H and K, bandpasses, respectively. A 
detailed description of the observing setup and the weather con- 
ditions can be found in Appendix A. 


! Database entry available at: 
cgi-bin/Dic-Simbad?/18721212 


http://cds.u-strasbg. fr/ 


Alexander J. Bohn et al.: Discovery of a directly imaged planet to the young solar analog YSES 2 


The data reduction was performed with PynPoint (version 
0.8.1; Stolker et al. 2019) and included basic processing steps 
such as dark and flat calibration, bad pixel cleaning, sky sub- 
traction, and correction for the instrumental distortion along the 
vertical axis of the detector. To remove the stellar halo that is 
affecting approximately the innermost 1"2 around the coron- 
agraph, we utilized an approach based on reference star dif- 
ferential imaging (RDI; Smith & Terrile 1984) in combination 
with principal component analysis (PCA; Amara & Quanz 2012; 
Soummer et al. 2012). As the parallactic rotation of our YSES 
observations is usually less than a few degrees, classical PSF 
subtraction schemes, such as angular differential imaging (ADI; 
Marois et al. 2006), perform much worse compared to this ref- 
erence library approach. This method of combined RDI plus 
PCA was already successfully employed to recover circumstel- 
lar disks in archival HST data (e.g., Choquet et al. 2014) and 
within the scope of our survey for the discovery of a transition 
disk around the YSES target Wray 15-788 (Bohn et al. 2019). 

Owing to the same location on sky, similar distances, and 
spectral types, all YSES targets exhibit very similar magnitudes 
in the red part of the optical spectrum (where the wavefront 
sensor of SPHERE is operating) and at the near-infrared wave- 
lengths of our scientific observations. This facilitates compara- 
ble AO corrections amongst all our YSES observations, and the 
resulting images compose an excellent reference library to per- 
form RDI. The reference targets that were used for our library 
PSF subtraction are listed in Appendix B. We modeled the stel- 
lar PSF with 50 principal components that were obtained from 
our full reference library. After the PSF subtraction, the frames 
were de-rotated according to their parallactic angles and median 
combined. 

For the astrometric calibration we used the standard in- 
strumental solution as presented by Maire et al. (2016) with 
a wavelength independent true north offset of —1°75 + 0°08 
and plate scales of (12.251 + 0.010) maspx^! and (12.265 + 
0.010) mas px^! in the H and K, band, respectively. 


4. Stellar properties 


We briefly summarize previous literature characterizing YSES 2 
in $4.1 and compile the stellar properties of YSES 2 in Table 
1. In $4.2 and $4.3, we derive updated stellar parameters for 
YSES 2, more importantly, including stellar mass and age. 


4.1. Previous studies 


YSES 2 was first identified as a young star in the Search for As- 
sociations Containing Young stars (SACY) survey (Torres et al. 
2006) of optical counterparts to the ROSAT All-Sky Survey 
(RASS) X-ray sources (Voges et al. 1999). Torres et al. (2006) 
reported the star to be a Li-rich (EW[LiI 46707] = 367 mÀ) 
K1V(e) star with filled-in Ha, showing fast rotation (vsini = 
19.3 km s^!) and radial velocity 15.8 + 1.0 km s™!. Based on its 
position, proper motion, and youth indicators, Preibisch & Ma- 
majek (2008) included the star in a list of new members of the 
LCC subgroup of the Sco-Cen OB association (their Table 4), 
and provided initial estimates of isochronal age (16 Myr), mass 
(1.1 Mo), and fractional X-ray luminosity (log(Lx/Lpo1) = -3.2). 
Preibisch & Mamajek (2008) also predicted a kinematic distance 
of 109 pc (based on the proper motion and space velocity of 
LCC), which compares remarkably well to the Gaia EDR3 paral- 
lactic distance (w = 9.1537 + 0.0118 mas, D = 109.25 + 0.14 pc; 
Gaia Collaboration et al. 2020). Kiraga (2012) reported the star 


to be a variable in the All Sky Automatic Survey (ASAS; IDed as 
ASAS J112755-6625.9) showing high amplitude (0.093 mag in 
V) and rapid rotation (Prot = 2.7325 d), which is consistent with 
the observed saturated X-ray emission (and right near the median 
rotation period for Sco-Cen pre-main-sequence stars of (Prot) ~ 
2.4 d; Mellon et al. 2017). Pecaut & Mamajek (2016) include 
the star in their age analysis of pre-main-sequence K stars across 
Sco-Cen, estimating an age of 23 Myr and a mass of 1.0 Mo. 
The star has subsequently appeared in multiple LCC member- 
ship lists (Gagné et al. 2018; Goldman et al. 2018; Damiani et al. 
2019). The status of this star as a pre-main-sequence member of 
LCC is strongly corroborated by a 99.9 % membership probabil- 
ity from the BANYAN X algorithm (Gagné et al. 2018) applied 
to the available Gaia astrometry and radial velocities (Gaia Col- 
laboration et al. 2018, 2020). 


4.2. Spectral analysis of YSES 2 


To check the previously published spectral properties of YSES 2, 
we examined two archival UVES spectra from the ESO archive 
taken on UT 2007 May 2 (Program 079.C-0556(A); PI Torres). 
The UVES spectra at resolution R = 40,000 were convolved to 
lower resolution R = 3,000 and compared to the grid of MK spec- 
tral standards from Pecaut & Mamajek (2016). The blue spec- 
trum, 3280 Å- 4560 A, is consistent with KOV, but the red spec- 
trum, 4730 Å- 6840 A, appears to be K2V. The Ha line exhibits 
marginal emission, similar to the filled-in emission reported by 
Torres et al. (2006). Hence, we infer a temperature type of K1+1 
and confirm the type K1V(e) published by Torres et al. (2006). 
From the original spectra, we independently measure the equiv- 
alent width of the Lil 46707 feature to be 364 +5 mA, by simul- 
taneously fitting Voigt profiles to the Lil feature and Fe I blend 
nearby (see, e.g., Soderblom et al. 1993). This is in good agree- 
ment with the 367 mÀ reported by Torres et al. (2006). The as- 
signed luminosity class that we adapted from Torres et al. (2006) 
does not necessarily imply that YSES 2 is on the main sequence 
rather than being a pre-main-sequence star. Even though the lu- 
minosity and gravity indicators used by Torres et al. (2006) were 
more in line with main-sequence dwarfs than subgiant or giant 
standard stars, more persuasive indicators such as the HRD po- 
sition, Li absorption, X-ray emission, and the confirmed LCC 
membership clearly favor the pre-main-sequence evolutionary 
stage. 


4.3. Updated stellar parameters 


Using the VOSA spectral energy distribution (SED) analyzer 
(Bayo et al. 2008), we fit synthetic stellar spectra to the ob- 
served visible and infrared photometry for YSES 2. For priors, 
we constrained the reddening to be E(B — V) = 0.016 x 0.017 
mag based on the STILISM 3D reddening maps (Lallement et al. 
2019) and searched for best-fit synthetic spectra in the range 
3000 K < Teg < 6000 K, 3.5 dex < log g < 4.5 dex, and metallici- 
ties -0.5 < [Fe/H] < 0.5 and [a/Fe] = 0.0. A Bayesian fit using the 
BT-Settl-CIFIST models using 22 photometric points yielded the 
following parameters: Ay = 0.06 (0.02-0.09; 68%CL; 0.00-0.11; 
95%CL), Te = 4749 K (4709-4789 K; 68%CL; 4700-4900 K; 
95%CL), log g = 3.9 dex (3.5-4.5 dex), [Fe/H] = 0.0. We present 
the results of this SED fit in Fig. 1. 

The best-fit bolometric flux is fj; = (1.7490 + 0.0248) x 107? 
ergs”! cm, which on the IAU 2015 apparent bolometric mag- 


? Online available at: http://svo2.cab.inta-csic.es/theory/ 
vosa/index.php 


Article number, page 3 of 15 


A&A proofs: manuscript no. ajbohn yses companion tyc8984 


Table 1. Stellar properties of YSES 2. 


Parameter Value Ref. 
Main identifier YSES 2 (1) 
TYCHO ID TYC 8984-2245-1 (2) 
2MASS ID J11275535-6626046 (3) 
GAIA EDR3 ID 523679288033301 1968 (4) 
a (J2000) [hh mm ss.sss] 11 27 55.355 (4) 
ó (J2000) [dd mm ss.ss] -66 26 04.50 (4) 
Spectral Type K1V(e) (5) 
w [mas] 9.1537 € 0.0118 (4) 
D [pc] 109.25 +0.14 (7) 
Hox [mas yr! ] -34.025 + 0.013 (4) 
us [mas yr! ] 2.319 € 0.011 (4) 
Vraa [km s^] 13.41 x 0.17 (5) 
B [mag] 11.819 + 0.010 (8) 
V [mag] 10.860 + 0.017 (8) 
G [mag] 10.525 + 0.003 (4) 
I [mag] 9.773 x 0.044 (9) 
J [mag] 9.006 + 0.026 (3) 
H [mag] 8.484 + 0.029 (3) 
K; [mag] 8.358 + 0.029 (3) 
W1 [mag] 8.323 + 0.014 (10) 
W2 [mag] 8.351 + 0.008 (10) 
W3 [mag] 8.258 + 0.019 (11) 
W4 [mag] 7.929 € 0.118 (11) 
Prot [day] 2.7325 (9) 
vsin i [km s^!] 19.3 € 0.5 (5) 
log(Lx/Lyo1) [dex] -3.07 + 0.23 (9) 
EW(Li 146707) [mA] 364 + 0.05 (6) 
EW(Ho) [mA] 0.0 (5) 
U [km s^! ] -10.10 + 0.08 (7) 
V [km s! -18.93 € 0.12 (7) 
W [km s7!] -5.60 + 0.09 (7) 
Ay [mag] 0.067007 (6) 
Ter [K] 4749 + 40 (6) 
my [mag] 10.396 + 0.015 (6) 
Mboi [mag] 5.204 + 0.016 (6) 
log (L/Lo) [dex] -0.1854 + 0.0063 (6) 
R [Ro] 1.193 + 0.022 (6) 
Mass [Mo] 1.10 x 0.03 (6) 
Age [Myr] 13.9 2.3 (6) 


References. (1) This paper, see $2; (2) Hgg et al. (2000); (3) Cutri et al. 
(2003); (4) Gaia Collaboration et al. (2020), and Gaia EDR3 and DR2 
ID £s are the same; (5) Torres et al. (2006); (6) this paper, see 84.3; 
(7) distance and heliocentric Galactic Cartesian velocity calculated us- 
ing Gaia EDR3 values (D = 1/@); (8) Henden et al. (2016); (9) Kiraga 
(2012); (10) Eisenhardt et al. (2020); (11) Cutri & et al. (2014). 


nitude scale translates to mp = 10.396 + 0.015 mag. Adopting 
the Gaia EDR3 parallax, this translates to absolute bolomet- 
ric magnitude Mgo = 5.204 + 0.016 mag and bolometric lumi- 
nosity log(L/Lo) = -0.1854 + 0.0063. This is considerably more 
accurate than previous estimates (log(L/Lo) = -0.06 (Preibisch 
& Mamajek 2008), log(L/ Lo) = -0.265 + 0.075 (Pecaut & Ma- 
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Fig. 1. Spectral energy distribution of YSES 2. The colored markers 
indicate the archival photometric measurements of the star and the blue 
curve presents our best-fit BT-Settl-CIFIST model to the data. The un- 
certainties of the photometric measurements are too small to be visual- 
ized in the figure. 


majek 2016)) and benefits from a very precise distance, well- 
constrained extinction from 3D reddening maps, and integrat- 
ing synthetic SEDs using 22 photometric data points. Combining 
this improved luminosity estimate with the improved Teg from 
the SED fitting (Teg = 4749+ 40K) yields a good estimate of 
the radius of the star (1.193 + 0.022 Ro). A comparison against 
evolutionary tracks from Baraffe et al. (2015) provided an up- 
dated stellar mass of (1.10 + 0.03) Mo and an age of (13.9 + 
2.3) Myr. Furthermore we note that the SED of the star showed 
no signs of infrared excess through to the WISE-4 band (22 um). 
The 2MASS-WISE colors (Ks-W1 = -0.035 + 0.032 mag, K,- 
W2 = -0.007 + 0.030 mag, K,-W3 = 0.100 + 0.035 mag, K;-W4 
= 0.429 + 0.122 mag) can be compared to the mean for K1 pre- 
main-sequence stars from Pecaut & Mamajek (2013) (K;-W1 
= 0.09 mag, K;-W2 = 0.06 mag, K;-W3 = 0.10 mag, K;,-W4 = 
0.18 mag), and show a significant hint of IR excess. The K;- 
W4 color is marginally red (2c excess), perhaps hinting at a 
debris disk (common among non-accreting pre-main-sequence 
Sco-Cen stars), but there is no corroborating evidence to further 
support this. 


We searched for common proper motion stellar or substel- 
lar companions for YSES 2. Assuming a mass of 1.0 Mo, the 
tidal radius of YSES 2 is ~1.35 pc (projected radius ~0°.71 or 
~2560”; Mamajek et al. 2013), that is, bound companions would 
be expected to lie projected within this radius. Surveying the lists 
of Sco-Cen candidates and pre-main-sequence stars from the 
Gaia DR2 catalogs of Goldman et al. (2018), Zari et al. (2018), 
and Damiani et al. (2019) with parallactic distances of «140 pc 
shows only a couple of likely LCC siblings within 1? of YSES 2: 
the poorly studied classical T Tauri star Wray 15-813 at 2845" 
(D = 101 pc) (Pereira et al. 2003) and uncharacterized candi- 
date pre-main-sequence object 2MASS J11375287-6631197 (D 
= 104 pc). A query of the recently released Gaia EDR3 catalog 
searching for co-moving, co-distant objects (with generous se- 
lection range of proper motions in œ and 6 within +5 mas yr !, 
and parallax +2 mas of YSES 2) yields zero candidate compan- 
ions within 1?. Thus far, YSES 2 appears to be a stellar singleton. 


Alexander J. Bohn et al.: Discovery of a directly imaged planet to the young solar analog YSES 2 


5. Observational results and analysis 


In $5.1 we show that our observations reveal a co-moving com- 
panion to YSES 2. The reduced images for both epochs are pre- 
sented in Fig. 2. Our photometric analysis in $5.2 indicates that 
this companion has a mass that is significantly lower than the 
deuterium burning limit of ~ 13 Mjyp. We refer to this newly 
identified planet as YSES 2b henceforth. In $5.3 we present the 
detection and mass limits of our acquired data. 


5.1. Companion astrometry 


We extracted the companion astrometry and photometry by the 
injection of negative artificial companions (e.g., Lagrange et al. 
2010; Bonnefoy et al. 2011). A detailed description of our 
method is presented in Appendix C. The extracted astrometry 
is listed in Table 2. As visualized in the proper motion plot 
in Fig. 3 the companion is clearly incompatible with the cal- 
culated trajectory of a static background object at 140 signifi- 
cance. The relative astrometric motion with respect to the pri- 
mary is consistent with a comoving companion. This conclusion 
is further confirmed by a similar analysis of other point sources 
within the detector field of view. As presented in Appendix D, 
all additional off-axis point sources are consistent with being 
non-moving background contaminants. We thus conclude that 
YSES 2b is a gravitationally bound companion to its solar-mass 
host star. From our astrometric measurements we derived a pro- 
jected physical separation of approximately 115 au. Future as- 
trometric measurements are required to derive meaningful con- 
straints to the orbital parameters of this wide-orbit planet. 


5.2. Companion photometry 


We present the photometry of the companion in Fig. 4 in a 
color-magnitude diagram. The corresponding numerical values 
are reported in Table 2. YSES 2b is consistent with a late L 
to early T spectral type when comparing it to colors of field 
brown dwarfs from the NIRSPEC Brown Dwarf Spectroscopic 
Survey (McLean et al. 2003, 2007), the IRTF spectral library 
(Rayner et al. 2009; Cushing et al. 2005), the L and T dwarf data 
archive Knapp et al. (2004); Golimowski et al. (2004); Chiu et al. 
(2006), and the SpeX Prism Libraries (Burgasser et al. 2010; 
Gelino & Burgasser 2010; Burgasser 2007; Siegler et al. 2007; 
Reid et al. 2006; Kirkpatrick et al. 2006; Cruz et al. 2004; Bur- 
gasser & McElwain 2006; McElwain & Burgasser 2006; Shep- 
pard & Cushing 2009; Looper et al. 2007; Burgasser et al. 2008; 
Looper et al. 2010; Muench et al. 2007; Dhital et al. 2011; Kirk- 
patrick et al. 2010; Burgasser et al. 2004). Object distances were 
derived from Gaia EDR3 (Gaia Collaboration et al. 2020), the 
Brown Dwarf Kinematics Project (Faherty et al. 2009), and the 
Pan-STARRSI 37 Survey (Best et al. 2018). In color-magnitude 
space, YSES 2b is very close to the innermost three planets of 
the HR 8799 multi-planetary system (Marois et al. 2008, 2010). 
These three planets are classified as mid to late L type dwarfs 
(e.g., Greenbaum et al. 2018), which agrees well with the se- 
quence evolution of the adjacent field brown dwarfs from L to 
T spectral types.? A similar spectral type in this domain, there- 
fore, seems very likely for YSES 2b, requiring confirmation by 


? We would like to note that the planets around HR 8799, although 
closely located to the sequence of field brown dwarfs in the selected 
SPHERE filters as presented in the color-magnitude diagram in Fig. 4, 
can have near-infrared colors in different passbands that are signifi- 
cantly distinct from those of their field dwarf analogs (e.g., Currie et al. 
2011). 


measurements at higher spectral resolution. Whereas the masses 
of the spectrally similar trio of HR 8799 c, d, and e are in the 
range 7-12 Mjyp (Wang et al. 2018; Marois et al. 2008, 2010), it is 
likely that YSES 2b has an even lower mass as the system age of 
(13.9+2.3) Myris significantly younger than the age of HR 8799, 
which is claimed to be member of the Columba association with 
an age of 30-50 Myr (Zuckerman et al. 2011; Bell et al. 2015). 
This is supported by the AMES-COND and AMES-dusty mod- 
els (Allard et al. 2001; Chabrier et al. 2000) that we present in 
Fig. 4 for a system age of 13.9 Myr. An individual evaluation of 
these isochrones yielded masses from 5.3 Mjup to 8.0 Mjyp as 
presented in Table 2. The uncertainties originate from the errors 
in the system age and planet magnitude that were propagated by 
à bootstrapping approach with 1,000 randomly drawn samples 
from Gaussian distributions around both parameters. When com- 
bining the posterior distributions for the different models and fil- 
ters we derived a final mass estimate of 6350 My as the 68% 
confidence interval around the median of the sample. This esti- 
mate is based on broadband photometric measurements alone; 
further spectral coverage of the planetary SED will be important 
to constrain its effective temperature, luminosity, surface gravity, 
and mass. 


5.3. Detection limits 


To derive upper mass limits for additional companions in the sys- 
tem, we calculated the detection limits of our datasets. As a base- 
line, we evaluated the contrast in the image that was obtained 
by de-rotating and median combining the individual exposures 
without any PSF subtraction. This image covers the full field of 
view of the SPHERE/IRDIS detector up to an angular separation 
of 5"5. We evaluated the contrast directly in the final image us- 
ing aperture photometry. The chosen aperture size was one full 
width at half maximum (FWHM) of the unsaturated stellar PSF 
as measured in the flux images (see Table A.1 in Appendix A). 
The signal flux was measured as the sum over the full circular 
aperture within the mean combined flux image and scaled for 
the flux difference with the science frames owing to the shorter 
exposure times and the applied neutral density filters. For sev- 
eral radial positions that were equidistantly sampled from 0/15 
to 57/50 in steps of 005, we measured the noise as the standard 
deviation of the integrated flux within apertures that were dis- 
tributed around the star at the same radial separation (excluding 
the signal aperture itself). We applied a sigma clipping with an 
upper bound of 3c to the integrated fluxes of the noise apertures 
before calculating the standard deviation to discard apertures that 
were polluted by flux of off-axis point sources (see the full frame 
image in Fig. D.1 in Appendix D). A correction for small sample 
statistics as described in Mawet et al. (2014) was considered in 
these noise calculations. We reiterated this analysis for six uni- 
formly spaced position angles and present the azimuthally aver- 
aged results as two-dimensional contrast curves in Fig. 5. The 
solid yellow lines represent the 5c raw contrast in the H and K; 
bands that was obtained without any PSF subtraction. 

For the innermost region around the star (< 172), the sensi- 
tivity was additionally assessed considering our PSF subtraction 
by RDI plus PCA. We used the ContrastCurveModule from 
PynPoint version 0.6.0 that utilizes the same aperture pho- 


^ As mentioned before, version 0.8.1 of PynPoint was used for all re- 
maining analysis steps. The modules of both versions are compatible; 
only the implementation of some algorithms changed throughout the 
development process. This affects the ContrastCurveModule, which 
follows the iterative process described in this paragraph for release ver- 
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Fig. 2. Multi-epoch observations of YSES 2 and its planetary-mass companion. Final data products of the SPHERE observations collected in the H 
band (left panel) and K, band (right panel) are presented. For both filters, the stellar PSF is modeled by 50 principal components that were derived 
from a reference library of YSES targets. These PSF models were subtracted and the residuals rotated such that north points up and east toward the 
left. In the presented images the median of these de-rotated residuals is shown. For the K, band data, uncorrected residuals of a wind-driven halo 
are detected that extend from the northeast to southwest. The planet YSES 2b is highlighted by white arrows. The primary is located at the origin 
of the coordinate system and we artificially masked the inner region up to the radial extent of the coronagraphic mask of 100 mas. To assess the 
spatial extent of the instrumental PSF, the median combination of the non-coronagraphic flux images of the primary star are shown in the lower 
left of each panel. The intensity of each flux image is rescaled to match the maximum and minimum counts in the corresponding residual science 


image, and we display both images with the same spatial and color scales. 


Table 2. Astrometry, photometry, and derived masses of YSES 2b. 


Observation date Filter Separation Position angle AMag Maps Photometric mass 
COND Dusty 
(yyyy-mm-dd) (mas) ©) (mag) (mag) (Miup) (Mp) 
2018-04-30 H 1057 +3 205.3 + 0.2 10.37 + 0.11 13.57 + 0.11 5.3 2 0.5 8.0 + 0.7 
2020-12-08 Ks 1053 +5 205.2 + 0.2 9.55+0.11 12.71+0.11 6.1+0.6 6.4+ 0.6 


tometry framework and metric to evaluate the contrast for sev- 
eral positions that were distributed around the star in the resid- 
ual images. For each position, the module injects an artificial 
companion, whose detection significance is evaluated after the 
PSF subtraction with RDI combind with PCA. In this frame- 
work, the signal aperture is directly placed on top of the posi- 
tion at which the artificial companion has been injected, and the 
noise apertures are azimuthally distributed around the primary 
star as described before, yet excluding the signal aperture itself. 
The companion template was obtained as the median combina- 
tion of the non-coronagraphic flux images that was scaled for 
the difference in exposure time and the attenuation due to the 
applied neutral density filter. From an initial magnitude contrast 
of 8 mag, the flux of the injected companion was adjusted and 
the post-processing was performed iteratively until the artificial 
companion was retrieved at 5c detection significance in the final 
image product. These limiting magnitude contrasts were stored 


sion 0.6.0. We prefer this implementation over the solution presented 
in PynPoint version 0.8.1, which calculates one attenuation factor per 
position that is based on the single injection of an artificial companion 
with a user-defined signal-to-noise ratio. 
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for each of the injection positions. For the calculation of the final 
contrast curves with RDI plus PCA, we used a radial sampling 
in the range 0715-1720 with a spatial resolution of 0'/05, and six 
position angles that were equidistantly sampled in polar space. 
Again, the contrast as a function of radial separation was ob- 
tained by azimuthal averaging of the various position angles. We 
considered several numbers of principal components to model 
the stellar PSF as indicated by the sequentially colored, dashed 
lines in Fig. 5. The detectable planet masses that correspond to 
the calculated magnitude contrasts were derived by evaluation of 
AMES-COND models at the system age of 13.9 Myr (see right 
axes of Fig. 5). 


The contrast performance close the star improves for an in- 
creasing number of principal components. This differential gain 
in contrast ceases for ~50 subtracted components and the con- 
trast for 100 principal components does not change significantly 
compared to the curve generated for half as many components. 
This justifies our previously selected value of 50 principal com- 
ponents that were used for our PSF subtraction with RDI plus 
PCA. This amount of components is equivalent to 19 % and 30 % 
of the reference libraries in the H and K, bands, which are com- 
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Fig. 3. Proper motion plot for YSES 2b. The colored markers represent 
the relative astrometry with respect to the primary star measured for our 
two observational epochs. The blue trajectory indicates the simulated 
motion of a static background object at infinity and the white marker 
is the theoretical position of such an object at our second observational 
epoch. 
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Fig. 4. Color-magnitude diagram for YSES 2b. The colored markers 
show the sequence of field brown dwarfs with various spectral types 
from M to late T. The white markers represent known directly imaged 
companions that are usually younger than the presented field objects. 
YSES 2b is highlighted by the red star. We further show AMES-COND 
and AMES-dusty evolutionary models that were evaluated at a system 
age of 13.9 Myr (solid lines). The markers along the line indicate the 
equivalent object masses in Mjyp. 


posed of 269 and 164 individual frames, respectively. In the H 
band we observe a contrast improvement of more than two mag- 
nitudes at an angular separation of 0"2. This corresponds to an 
increase in planet detection sensitivity by more than 45 Mjyp at 
this close separation. The contrast improvement in the H band 
is maximized at an angular separation of ~0’’3, where RDI plus 
PCA provides detection limits that are approximately three mag- 
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Fig. 5. 5c detection limits of our SPHERE/IRDIS observations in the 
H (upper panel) and K, bands (lower panel). On the left axis the magni- 
tude contrast with respect to the primary star is reported, and the abso- 
lute magnitudes are converted to detectable planet masses with AMES- 
COND models as indicated on the right axis; this scale varies between 
the H and K; bands. The solid yellow lines represent the limits when 
no PSF subtraction is performed. The dashed lines indicate the sensitiv- 
ity, when a PSF subtraction with RDI plus PCA is performed. The red 
star highlights the contrast of YSES 2b that we detect at ~ So signifi- 
cance in both filters after PSF subtraction with more than 50 principal 
components. 


nitudes deeper than our raw data. At angular separations larger 
than 1" the contrast improvement decreases as the flux con- 
tribution of the stellar PSF becomes negligible. At separations 
> 2" we reach a fundamental noise floor that is mainly com- 
posed of residual sky background and detector read out noise. 
The K; band contrast behaves very similar to the detection lim- 
its in the H band and the RDI plus PCA reduction scheme can 
provide a maximum gain of up to 2.5 mag at an angular separa- 
tion of ~0”4. The overall improvement for separations < 1" is 
marginally worse compared to the H band data and the contri- 
bution of the asymmetric wind-driven halo is clearly visible for 
separations in the range 075-172. Combining the data from the 
H and K; bands allows us to exclude stellar and brown dwarf 
companions around YSES 2 with masses > 13 Mj, for angular 
separations that are larger than 015. At 0"5 we are sensitive to 
objects that are more massive than 6 Mj, and for angular sepa- 
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rations that are larger than 2" we can even rule out planets with 
masses as low as 2 Myyp . 

This demonstrates, impressively, how a large reference li- 
brary can help to significantly improve the contrast performance 
at small angular separations < 1"2. Especially for datasets with 
little parallactic rotation, RDI plus PCA should be considered as 
a default PSF subtraction strategy. This conclusion is also sup- 
ported by first results from the star-hopping mode that was re- 
cently implemented at VLT/SPHERE (Wahhaj et al. 2021). As 
visualized by the red stars in Fig. 5, RDI plus PCA is required 
to detect YSES 2b at 5c significance in both the H and K, band 
data. 


6. Discussion 


The newly discovered planetary companion to YSES 2 is among 
the lowest mass direct imaging companions known to date. The 
only objects of similar low mass are 51 Eri b (Macintosh et al. 
2015), HD 95086b (Rameau et al. 2013), HR8799b (Marois 
et al. 2008), PDS 70 b and c, (Keppler et al. 2018; Haffert et al. 
2019; Wang et al. 2020; Stolker et al. 2020a), and YSES 1c 
(Bohn et al. 2020b). Of these, only YSES 1c is located around 
a solar-mass star. Within the uncertainties, the mass of YSES 2b 
is the same as YSES 1c, which we previously discovered in our 
survey. Given our mass estimate for the planet, the mass ratio of 
YSES 2b to its host star is g = 054. %. This value is com- 
parable, but slightly lower than the g = 0.57 + 0.10 % derived 
for YSES 1c. The mass ratio is the lowest among direct imaging 
companions to solar-type stars.? 

The in situ formation of super-Jovian planets at tens or hun- 
dreds of Astronomical Units is challenging. We recently dis- 
cussed possible formation scenarios for such objects in the con- 
text of the YSES 1 system in Bohn et al. (2020a). In our pre- 
vious study we considered scattering or planet capturing events 
to explain the current large separation of YSES 1b. However, 
dynamic scattering by a third body in the system is expected 
to produce high eccentricities, inconsistent with orbital stabil- 
ity of both planetary components in the YSES 1 system (Bohn 
et al. 2020b).° The new detection of YSES 2b in our small survey 
sample of 70 solar-type systems in Sco-Cen makes the hypoth- 
esis that we see captured free-floating planets unlikely. Goulin- 
ski & Ribak (2018) find with numerical simulations that only 
~0.1% of solar-type stars in the Galactic thin disk should cap- 
ture a free-floating planet in their lifetime. YSES 2b in principle 
might have formed in situ via disk gravitational instability. Boss 
(2011) find that they can produce 1-5 Mjyp planets between 30 au 
and 70 au with eccentricities as high as 0.35. If YSES 2b is on 
an eccentric orbit it may explain its current projected separation 
of 110 au. Kratter et al. (2010) conversely find with their hy- 
drodynamic simulations that planets formed via disk instability 
need to be at large separations outside of the 40 au to 70 au range 
to not accumulate too much mass and remain in the planetary 
regime. This hypothesis that gravitationally instabilities predom- 
inantly create brown dwarf and stellar companions is supported 
by other theoretical studies (e.g., Zhu et al. 2012; Forgan & Rice 
2013). Spatially resolved observations of gas-rich planet form- 
ing disks in the last few years have shown that radial substruc- 
tures, which are thought to be caused by perturbing planets, are 


> See Bohn et al. (2020a) for an overview of mass ratios of direct imag- 
ing companions to solar-type stars. 

5 We, however, note that in some cases high eccentricities are inferred 
for substellar companions, for example, the brown dwarf companion 
to the young solar analog PZ Tel (Mugrauer et al. 2012; Ginski et al. 
2014). 
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nearly ubiquitous (see, e.g., Garufi et al. 2018 for an overview, 
observed in scattered light). These structures can, in a growing 
number of cases, be traced out to tens or hundreds of Astro- 
nomical Units (e.g., ALMA Partnership et al. 2015; Ginski et al. 
2016; de Boer et al. 2016; van Terwisga et al. 2018). Recently it 
was found that these substructures are already present in proto- 
stellar disks as young as ~0.1 Myr (Sheehan et al. 2020; Sheehan 
2020), suggesting that planet formation sets in early and operates 
on short timescales. Despite this abundance of substructures ob- 
served in young, circumstellar environments, it is unlikely that 
all these protoplanetary disks are gravitationally unstable and 
support planet formation via this channel (e.g., Kratter & Lodato 
2016). Conversely, we expect that the timescale to form a planet 
via core accretion at the current location of YSES 2b would be 
too long, given the system age and that the gas-rich disk in the 
system has already dissipated (e.g., Haisch et al. 2001). 

Even though capturing scenarios are considered unlikely for 
YSES 2b, we cannot confidently conclude whether its forma- 
tion via either top-down or bottom-up scenarios is more likely: 
whereas a mass of 6.3 Mj, is rather low for an object to orig- 
inate from gravitational disk instabilities, core accretion would 
favor a formation at closer separations to the star. More data are 
thus required to explore the origin of this wide-orbit Jovian giant. 
Promising methods to evaluate the likelihood of either forma- 
tion scenario are by characterization measurements of the plan- 
etary atmosphere, continuous orbital monitoring to constrain es- 
pecially its eccentricity, and deeper searches for additional com- 
panions in the system. 

To identify the formation channel of YSES 2b via atmo- 
spheric characterization, we can utilize the framework postulated 
by Óberg et al. (2011), who argued that elemental abundances in 
the planetary atmosphere (and especially the C/O ratio) are di- 
rectly linked to its natal environment in the planet-forming disk. 
Different ice lines in the protoplanetary disk and the associated 
freeze out of the corresponding molecular species alter molecu- 
lar abundance ratios as a function of radial separation from the 
host star, making this atmospheric quantity a promising indica- 
tor of the natal environment and formation channel of a planet. 
The chemical and dynamical evolution of the disk can alter these 
initial abundances and should be considered in the analysis (e.g., 
Ali-Dib et al. 2014; Mordasini et al. 2016; Eistrup et al. 2016, 
2018). Gravity Collaboration et al. (2020) utilized this frame- 
work to study £ Pic b and proposed its formation through core 
accretion, with strong planetesimal enrichment based on its sub- 
solar C/O abundance ratio. If YSES 2b formed via disk grav- 
itational instability, then we expect this object to have similar 
elemental abundances as the primary star in the system, while 
formation by core accretion should lead to an over-abundance 
of heavy elements due to pebble accretion. YSES 2b, along with 
other planet mass objects detected by direct imaging, provides an 
ideal test case for future detailed atmospheric characterization. 

We can further continue to monitor the separation and po- 
sition angle of YSES 2b with respect to the primary star to de- 
rive orbital solutions for the planet (e.g., Wang et al. 2018). In 
particular, the VLTI/GRAVITY instrument (Gravity Collabora- 
tion et al. 2017) will be extremely useful for this purpose, as it 
facilitates an unprecedented astrometric precision down to sub- 
milliarcsecond scales (e.g., Gravity Collaboration et al. 2019; 
Wang et al. 2021). These astrometry measurements could be 
complemented by VLT/CRIRES* data to constrain the radial ve- 
locity of the planet and to obtain three-dimensional information 
about its orbital motion (e.g., Schwarz et al. 2016). The eccen- 
tricity of the planet might provide hints regarding the likelihood 
of a potential migration of the companion, which would be an 
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indicator of formation via core accretion at closer separation to 
the star. If this migration was caused by scattering off another, 
so far undetected companion to the primary star, a deep imaging 
campaign is required to search for evidence of such an additional 
component to the planetary system. At the moment, we cannot 
provide conclusive evidence for the most likely formation sce- 
nario of YSES 2b based on the available data; but future obser- 
vations might be able to shed light on the origin of this Jovian 
gas giant. 

Even though the second epoch observations of YSES are 
not concluded yet and candidate companions to —45 stars of 
our sample need to be confirmed or rejected, our survey has al- 
ready discovered three planetary-mass companions amongst 70 
young, Sun-like stars. This high planet-detection rate is in stark 
contrast to previous surveys that were targeting Sun-like stars at 
closer distances than the LCC (e.g., Kasper et al. 2007; Biller 
et al. 2013; Galicher et al. 2016), which discovered mostly stel- 
lar and brown dwarf companions. These preliminary statistical 
results from YSES tentatively indicate that despite the farther 
distance, Sco-Cen and especially LCC are more favorable than 
moving groups in the immediate solar neighborhood for the de- 
tection of young planets briefly after their formation. Many of 
the moving groups that were targeted in these aforementioned 
surveys are significantly older than LCC (15 + 3 Myr), such as 
the Tucana-Horologium moving group (45 + 4 Myr Bell et al. 
2015), the AB Dor moving group (149*75 Myr; Bell et al. 2015), 
or the Hercules-Lyra association (257 + 46 Myr; Eisenbeiss et al. 
2013). As a consequence of the decreasing luminosity of objects 
below the deuterium burning limit with increasing age (e.g., Bur- 
rows et al. 1997), it is natural that the sensitivity to Jovian plan- 
ets is worse around members of these associations compared to 
significantly younger host stars. Yet some of these closer mov- 
ing groups have ages comparable to that of LCC — such as the 
TW Hya association (10 + 3 Myr) or the 8 Pic moving group 
(24 + 3 Myr; Bell et al. 2015) — and should provide even bet- 
ter planet-detection sensitivities owing to their much closer dis- 
tances. However, before speculating about potential reasons for 
this tentative overabundance of planetary-mass companions to 
our YSES targets, it is necessary to finish the second epoch ob- 
servations of the survey and to derive reliable occurrence rates 
of planetary-mass companions to Sun-like stars in Sco-Cen. 


7. Conclusions 


We report the detection of a new directly imaged planet to the 
solar-mass primary YSES 2 that was discovered within the scope 
of YSES. Reassessment of the stellar parameters provided an 
effective temperature of Teg = (4749 + 40) K, a luminosity of 
log (L/Lo) = —0.1854 + 0.0063, a mass of (1.10 + 0.03) Mo, and 
a system age of (13.9 + 2.3) Myr. 

We detect YSES 2b in two consecutive epochs collected on 
2018 April 30 and 2020 December 18 with VLT/SPHERE. The 
companion has a projected separation of approximately 1"/05, 
which translates to a physical minimum distance of ~115 au with 
respect to the primary star. Photometric measurements in the H 
and K; bands constrain a planet mass of du. Myyp according 
to AMES-COND and AMES-dusty evolutionary models. This 
mass estimate is supported by the position of the object in color- 
magnitude space, where it is located amongst the mid to late L 
type field brown dwarfs and close to HR 8799 c, d, and e. The 
slightly higher mass estimates of these exoplanets on the order of 
7—12 My, are consistent with the older system age of HR 8799 
of 30-50 Myr. 


The mass and separation of YSES 2b are inconsistent with 
planet populations for most in situ formation scenarios: whereas 
disk instabilities predominantly create companions above the 
deuterium burning limit at a separation of 110 au, core-accretion 
mechanisms are not efficient enough to form a planet of 6.3 Mjup 
this widely separated from the primary star. So, the new com- 
panion might be either at the low-mass end of potential in situ 
formation outcomes from top-down scenarios, or it formed via 
core accretion at closer separation to the star and migrated to its 
current location. Atmospheric characterization measurements of 
molecular abundance ratios, orbital monitoring, and evaluation 
of the eccentricity of the planet, or a deep search for additional 
companions in the system, might help to evaluate the likelihood 
of these potential formation pathways. While we cannot rule out 
scattering or capture scenarios, we point out that the former re- 
quire an (as of yet) third undetected body in the system, while 
the latter are unlikely given numerical simulations. YSES 2b is 
an important addition to the sparsely populated group of wide- 
orbit gas giant companions. Owing to the moderate separation 
with respect to the primary star, spectroscopic observations with 
JWST, VLT/ERIS, or VLTI/GRAVITY will be easily available. 
These data will be important to further constrain the properties of 
this Jovian companion. Measurements of molecular abundance 
ratios such as C/O or its orbital eccentricity might even facilitate 
hypotheses regarding the most likely formation mechanism for 
this wide-orbit gas giant planet. 


Our data rule out brown dwarf and stellar companions with 
M > 13 Mjyp in the SPHERE/IRDIS field of view for angular 
separations > 0715 and at 0" 5 we can exclude objects that are 
more massive than 6 Mj, . At separations that are larger than 
2" we are even sensitive to planets with masses as low as 2 Mjup 
. In general, the applied PSF subtraction scheme based on RDI 
plus PCA is extremely successful and provides substantial con- 
trast improvements (> 1 mag) for separations that are smaller 
than 1". In the H band, the PSF subtraction enhances our sensi- 
tivity by more than 45 Mj, at 0"2, and the greatest contrast im- 
provement of ~3 mag is achieved at an angular separation of 0/3. 
Our YSES strategy with short snapshot observations of < 5 min 
combined with a large reference library for PSF subtraction is 
certainly a promising approach to image planetary-mass com- 
panions to young, Sun-like stars in Sco-Cen. With three newly 
discovered planetary-mass companions in less than 40h of al- 
located telescope time the survey efficiency is unprecedented 
and the mission concept can certainly be applied to future high- 
contrast imaging studies targeting different samples of pre-main- 
sequence stars. 
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Appendix A: Observational conditions and setup 


We present the observational setup and the weather conditions 
for our SPHERE observations in Table A.1. 


Appendix B: Reference library 


The reference libraries were compiled from the full amount of 
YSES data that were collected under ESO IDs 099.C-0698(A) 
(PI: Kenworthy), 0101.C-0153(A) (PI: Kenworthy), 0101.C- 
0341(A) (PI: Bohn), and 106.20X2.001 (PI: Vogt). We used RDI 
in the innermost region of the images < 1"2, where the stellar 
halo was dominating the received flux. We deselected all targets 
with obvious point sources or extended structures in this region 
because these signals are not part of the stellar PSF and would 
therefore deteriorate the quality of our model created by PCA. 
The remaining targets, their observation epochs, and observing 
conditions that were used as a reference library for the H and K, 
band data are listed in Tables B.1 and B 2, respectively. In the 
H band we have 269 individual reference frames and in the K; 
band we have 164. 


Appendix C: Extraction of companion astrometry 
and photometry 


We extracted the companion astrometry and photometry with the 
SimplexMinimizationModule of PynPoint. This injects an 
artificial planet into the data prior to the stellar PSF subtraction 
with RDI plus PCA. The planet template PSF is obtained from 
the unsaturated, non-coronagraphic flux images that were taken 
alongside the observations. The methods injects the artificial 
planet into the data at the approximate position and magnitude of 
the real point source, considering the parallactic rotation during 
the observing sequence. The PSF subtraction is performed using 
the same library as before (see Appendix B), we smooth the im- 
age with a Gaussian kernel with a FWHM of 12 mas (which cor- 
responds to the size of a detector pixel) to reduce pixel-to-pixel 
variations, and we evaluate the residuals in an aperture with a 
diameter of ~ 07/25 around the injection position. We choose the 
image curvature, which is represented by the determinant of the 
Hessian matrix as function of merit, which we aim to minimize 
by varying the input separation, position angle, and magnitude 
contrast of our artificial companion. We do not use the absolute 
value norm as presented by Wertz et al. (2017) as an objective 
to the minimization because this would not consider large-scale 
features in the residual image that are not correctly modeled by 
our PSF subtraction approach. Such a feature is for instance the 
asymmetric wind driven halo (Cantalloube et al. 2018) that is 
apparent in the K; band data in the northeastern to southwestern 
direction (see right panel of Fig. 2). This uncorrected stellar flux 
contributes to the planet signal and minimization of the absolute 
value norm around the planet position would certainly overesti- 
mate its flux and perhaps even compromise its astrometry. Planet 
separation, position angle, and magnitude contrast are optimized 
simultaneously by a Nelder-Mead simplex minimization algo- 
rithm (Nelder & Mead 1965). 

Owing to the optimization process our final values for the 
planetary astrometry and photometry do not exhibit any intrin- 
sic uncertainties. To derive the systematic uncertainties of our 
injection and minimization approach, we follow the analysis de- 
scribed by Stolker et al. (2020b), using the cube in which the op- 
timized negative planet is injected such that no companion signal 
remains in the data. For 24 position angles that are equidistantly 
distributed in polar space we inject positive artificial companions 
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Fig. D.1. Reduced SPHERE data for the full IRDIS field of view. The 
images are de-rotated and median combined; an unsharp mask is ap- 
plied to remove the stellar halo. Four additional companion candidates 
to YSES 2b are identified in the field of view. The image is presented 
at an arbitrary logarithmic color scale to highlight the off-axis point 
sources. 


into the data using the same magnitude contrast and the same ra- 
dial separation as previously determined for our companion. We 
extract the astrometry and photometry of these artificial compan- 
ions with the same method as described before and we evaluate 
the deviations from the injection position and flux. The standard 
deviation along the 24 distinct positions is utilized as uncertainty 
of our extraction method. These are combined with additional as- 
trometric uncertainties originating from the detector plate scale, 
the true north offset, and the centering accuracy of 2.5 mas (see 
SPHERE manual) to derive the final value of planet separation 
and position angle as presented in Table 2. For the companion 
photometry, we add uncertainties due to the variation of the un- 
saturated stellar PSF throughout the sequence of flux measure- 
ments and we account for transmissivity variations of the neutral 
density filter across the broadband filter that was used for our 
Observations (either the H or K; band). 


Appendix D: Astrometric analysis of background 
objects 


In addition to YSES 2b, there are four candidate companions 
(CCs) in the SPHERE/IRDIS field of view that we could iden- 
tify in both observational epochs. These CCs are presented in 
Fig. D.1, in which we show the de-rotated data from the night of 
2020 December 12. No PSF subtraction with RDI is performed, 
instead we just applied an unsharp mask with a Gaussian kernel 
size of 5 pixels. YSES 2b can easily be identified in this im- 
age product as well. For the remaining CCs, we present the rela- 
tive astrometric offsets between both observational epochs in the 
proper motion diagram in Fig. D.2. The relative motions of CCs 
2-5 are clearly compatible with stationary background objects, 
and co-movement can be ruled out for all of them. 
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Table A.1. SPHERE observations of YSES 2. 


Observation date Filter FWHM’ NEXPXNDITxDIT^ An‘ mu (XY (tof 
(yyyy-mm-dd) (mas) (1x1 xs) ©) ) (ms) 
2018-04-30 H 50.5 4x1x32 0.98 0.87 1.343 6.25 
2020-12-08 Ks 61.7 1x18x16 1.37 0.55 1.52 4.30 


Notes. © Full width at half maximum measured for the non-coronagraphic stellar PSF. © NEXP describes the number of exposures, NDIT is 
the number of subintegrations per exposure, and DIT is the detector integration time of an individual subintegration. © Az describes the amount 
of field rotation during the observation, if it is carried out in pupil-stabilized mode (only valid for CI observations). (^ (X) denotes the average 
airmass during the observation. “ (w) denotes the average seeing conditions during the observation. (ro) denotes the average coherence time 


during the observation. 
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Fig. D.2. Proper motion plot for background objects in the 
SPHERE/IRDIS field of view. The pink markers indicate the relative 
astrometric offsets to the first observational epoch that is plotted at the 
origin of the coordinate system (orange marker). The blue trajectory 
represents the simulated motion of a static background object at infin- 
ity and the white marker shows the relative positional offset of such an 
object at the time of our second observation. 
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Table B.1. Reference library for the data reduction in H band. 


Target Observation date NEXPXNDITxDIT^ (wy? (X) (19)d 
(2MASS ID) (yyyy-mm-dd) (1x1xs) (^) (ms) 
J11272881-3952572 2017-04-18 4x1x32 1.51 1.10 1.40 
J11320835-5803199 2017-06-17 4x1x32 0.67 1.47 2.90 
J11445217-6438548 2018-05-14 4x1x32 0.72 1.31 2.38 
J11454278-5739285 2018-06-04 4x1x32 0.70 1.19 2.80 
J11454278-5739285 2019-01-13 4x1x32 1.14 1.62 3.83 
J12065276-5044463 2017-04-02 3x1x32 1.24 1.11 1.50 
J12090225-5120410 2018-05-15 4x1x32 0.86 1.12 2.70 
J12090225-5120410 2019-12-14 12x2x32 0.63 1.51 7.15 
J12101065-4855476 2017-04-18 4x1x32 1.71 1.15 1.40 
J12113142-5816533 2018-12-22 3x2x32 1.46 1.47 2.13 
J12113142-5816533 2019-02-18 4x2x32 0.45 1.23 14.30 
J12160114-5614068 2018-12-27 4x2x32 0.41 1.45 11.88 
J12164023-7007361 2018-12-23 3x1x32 0.98 1.59 2.93 
J12164023-7007361 2019-02-15 4x1x32 0.54 1.63 11.20 
J12185802-5737191 2017-06-17 2x1x32 0.72 1.22 2.70 
J12195938-5018404 2018-12-30 4x1x32 0.53 1.62 8.00 
J12210499-7116493 2019-01-12 4x2x32 0.80 1.53 4.25 
J12220430-4841248 2017-04-18 3x1x32 1.82 1.17 1.40 
J12234012-5616325 2017-06-17 4x1x32 0.62 1.72 3.45 
J12264842-5215070 2018-12-30 4x1x32 0.40 1.38 8.20 
J12302957-5222269 2018-12-30 4x1x32 0.38 1.33 9.85 
J12333381-5714066 2019-01-01 4x1x32 0.76 1.37 7.03 
J12333381-5714066 2019-01-14 4x1x32 1.26 1.21 2.45 
J12361767-5042421 2018-12-30 4x1x32 0.51 1.59 4.68 
J12361767-5042421 2019-12-18 16x2x32 1.14 1.57 3.07 
J12374883-5209463 2018-12-30 4x1x32 0.41 1.50 7.30 
J12383556-5916438 2019-01-03 4x1x32 0.52 1.59 13.90 
J12383556-5916438 2019-01-12 4x1x32 0.79 1.26 4.25 
J12393796-5731406 2017-06-17 4x1x32 0.64 1.77 3.83 
J12404664-5211046 2018-04-30 4x1x32 0.74 1.13 7.05 
J12442412-5855216 2017-06-17 4x3x32 0.71 1.37 2.67 
J12454884-5410583 2018-04-30 4x1x32 0.71 1.15 6.92 
J12480778-4439167 2017-06-17 4x2x32 0.90 1.34 2.75 
J12505143-5156353 2019-01-12 4x1x32 1.14 1.32 3.75 
J12510556-5253121 2019-01-08 4x1x32 0.58 1.68 3.90 
J13015069-5304581 2019-01-08 4x1x32 0.55 1.60 3.95 
J13055087-5304181 2018-06-11 4x1x32 0.82 1.14 1.95 
J13055087-5304181 2018-07-04 4x1x32 1.73 1.14 1.70 
J13064012-5159386 2018-04-30 4x1x32 0.56 1.13 8.15 
J13065439-4541313 2018-04-08 4x1x32 0.46 1.09 5.65 
J13095880-4527388 2018-05-01 4x1x32 1.08 1.07 2.70 
J13103245-4817036 2018-05-01 4x1x32 1.03 1.09 3.30 
J13121764-5508258 2017-08-31 4x1x32 0.68 2.22 4.42 
J13121764-5508258 2018-05-15 4x1x32 0.62 1.16 2.50 
J13174687-4456534 2018-05-28 4x1x32 0.70 1.07 4.33 
J13334410-6359345 2017-07-05 4x1x32 1.06 1.53 3.05 
J13343188-4209305 2017-04-02 4x1x32 1.14 1.21 1.70 
J13354082-4818124 2017-04-02 4x1x32 1.06 1.30 2.08 
J13380596-4344564 2017-04-02 4x1x32 1.05 1.32 2.40 
J13455599-5222255 2018-04-28 4x1x32 0.64 1.13 6.35 


Notes. “ NEXP describes the number of exposures, NDIT is the number of subintegrations per exposure, and DIT is the detector integration time 
of an individual subintegration.  (X) denotes the average airmass during the observation. © (w) denotes the average seeing conditions during 
the observation. (? (ro) denotes the average coherence time during the observation. 
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Table B.2. Reference library for the data reduction in K, band. 


Target Observation date NEXPXNDITXDIT? (wy? (X)* (To)! 
(2MASS ID) (yyyy-mm-dd) (1x1xs) (^) (ms) 
J11445217-6438548 2018-05-14 4x1x32 0.77 1.31 2.60 
J11454278-5739285 2019-01-13 4x1x32 1.18 1.59 3.58 
J12090225-5120410 2018-05-15 4x1x32 0.70 1.12 2.90 
J12113142-5816533 2018-12-22 4x2x32 1.38 1.44 2.05 
J12113142-5816533 2019-02-18 4x2x32 0.45 1.22 15.00 
J12160114-5614068 2018-12-27 4x2x32 0.47 1.42 10.27 
J12164023-7007361 2018-12-23 4x1x32 1.06 1.58 3.43 
J12164023-7007361 2019-02-15 4x1x32 0.57 1.61 10.75 
J12195938-5018404 2018-12-30 4x1x32 0.55 1.59 9.00 
J12210499-7116493 2019-01-12 4x2x32 0.82 1.52 4.40 
J12264842-5215070 2018-12-30 4x1x32 0.41 1.36 9.20 
J12302957-5222269 2018-12-30 4x1x32 0.45 1.32 7.48 
J12333381-5714066 2019-01-01 4x1x32 0.80 1.36 6.25 
J12333381-5714066 2019-01-14 4x1x32 1.24 1.21 2.30 
J12333381-5714066 2020-12-10 1x20x16 0.58 1.76 5.50 
J12361767-5042421 2018-12-30 4x1x32 0.47 1.56 6.22 
J12374883-5209463 2018-12-30 4x1x32 0.46 1.48 6.95 
J12383556-5916438 2019-01-03 4x1x32 0.46 1.56 12.47 
J12383556-5916438 2019-01-12 4x1x32 0.94 1.26 3.45 
J12404664-5211046 2018-04-30 4x1x32 0.87 1.13 7.10 
J12454884-5410583 2018-04-30 4x1x32 0.66 1.15 8.97 
J12505 143-5 156353 2019-01-12 4x1x32 1.03 1.31 4.10 
J12510556-5253121 2019-01-08 4x1x32 0.52 1.65 3.98 
J13015069-5304581 2019-01-08 4x1x32 0.49 1.58 4.80 
J13055087-5304181 2018-06-11 4x1x32 0.93 1.14 2.02 
J13055087-5304181 2018-07-04 4x1x32 1.73 1.14 1.70 
J13064012-5159386 2018-04-30 4x1x32 0.56 1.13 9.88 
J13065439-4541313 2018-04-08 4x1x32 0.55 1.09 4.68 
J13095880-4527388 2018-05-01 4x1x32 1.03 1.07 2.45 
J13103245-4817036 2018-05-01 4x1x32 0.87 1.10 4.40 
J13121764-5508258 2018-05-15 4x1x32 0.62 1.16 3.00 
J13174687-4456534 2018-05-28 4x1x32 0.67 1.07 4.15 
J13455599-5222255 2018-04-28 4x1x32 0.65 1.13 6.03 


Notes. “ NEXP describes the number of exposures, NDIT is the number of subintegrations per exposure, and DIT is the detector integration time 
of an individual subintegration. ( (X) denotes the average airmass during the observation. (^ (w) denotes the average seeing conditions during 
the observation. (? (ro) denotes the average coherence time during the observation. 
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